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Our main concern in this work is to show how higher shell admixtures affect the spectrum of a 
Q ■ Q interaction. We first review how, in the valence space, the familiar SU(3) result for the energy 
spectrum can be obtained using a coordinate space Q Q interaction rather than the Elliott one which 
is symmetric in r and p. We then reemphasize that the Elliott spectrum goes as L(L + 1) where L is 
the orbital angular momentum. While in many cases this is compatible with the rotational formula 
which involves 1(1 + 1), where I is the total angular momentum, there are cases, e.g. odd-odd nuclei, 
where there is disagreement. Finally, we consider higher shell admixtures and devise a scheme so as 
to obtain results, with the Q ■ Q interaction, which converge as the model spaces are increased. We 
consider not only ground state rotational bands but also those that involve intruder states. 



On 



I. INTRODUCTION 



There have been studies in the past of the quadrupole-quadrupole (Q ■ Q) interaction in multishell spaces using 
the non-compact symplectic groups [ffl,pj . Most of the focus of these studies has been on the ground state rotational 
bands, the 7 vibration bands (K = 2), the giant quadrupole resonances and the effective charges. In this work we 
wish to discuss the problems of extending the above studies to other topics such as intruder states, odd-odd nuclei, 
etc. We feel that there has not been any clear discussion of how higher shell admixtures affect the overall low lying 
spectrum in a given nucleus, nor of how the low lying bands of odd-odd nuclei are described with this interaction. 
Q\ . Also we have recently shown how to restore Elliott's SU(3) results || with a coordinate space quadrupole-quadrupole 
0^ ' interaction rather than one that is symmetric in position and momentum. We will show that this is of more than 
academic interest. We also discuss differences between Elliott's SU(3) approach [|[j6) and the rotational model. The 
+1^ 1 problem of convergence with a Q ■ Q interaction is discussed in the last section. 

1— 1 I The Q ■ Q interaction has a long history, the early part of which can be found by reading excerpts from the text 
book by A.M. Lane 0, and for recent applications the book by I. Talmi is recommended ||. 

Following the notation of Golin and Zamick j^] we define the multipole-multipole interaction (including isospin) as 

V(r) = - X [(2L + 1) (2T + 1)] 1/2 [Y L (l) Y L (2)] ° / ( n ) / (r 2 ) \s TA [r (1) r (2)]° + S T . ] (I) 

!><!■- 

where -y/3 [t (1) r (2)]° = r (1) • r (2) . 

The detailed expressions for the unnormalized antisymmetrized matrix elements of this interaction are given in the 
above mentioned work of Golin and Zamick ||. The expression is fairly complicated, involving Racah coefficients. 
However the expression for the direct part of the particle-hole matrix element (ph) for a ph state with total angular 
momentum / and isospin T' is much simpler. The expression is proportional to 61 l6t,t' and indeed there is no other 
/, T' dependence. Thus, if one ignores the exchange part of the interaction one has a simple schematic interaction. 
Thus, as Lane describes Q, we can use this interaction to describe collective vibrational spectra. 

The pairing plus quadrupole model was very popular and has been used by Kisslinger and Sorensen |l0| , Baranger 
and Ikeda et al. around 1960, especially to describe the first excited 2 + state in vibrational nuclei. It can also 
be applied to giant quadrupole resonances as emphasized by Bohr and Mottelson Here, the schematic model of 
Brown and Bolsterli can be used. One can also couple a particle to the giant quadrupole resonance and get an 
E2 effective charge. If one uses the selfconsistent strength for \ as given by Bohr and Mottelson one finds that the 
isoscalar polarization charge Se p + 5e n is equal to the bare charge in an RPA calculation. 

In a different direction Elliott has shown || that a suitably defined Q ■ Q interaction leads to rotational spectra in an 
open shell nucleus with single-particle states defined by an isotropic three-dimensional harmonic oscillator potential. 
His interaction is 
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Note that it includes i = j terms and that his quadrupole operator is symmetric in position and momentum and this 
ensures that there will be no mixing of major shells. This will be discussed in more detail later. Elliott noted that 
the Casimir operator for the SU(3) group is 

C 2 = Q-Q-3L-L (3) 

From this he obtained the famous expression for the energy levels for a rotational band for states with SU(3) quantum 
numbers A, /i 

(Ve) = X [-4 (A 2 + + A^ + 3 (A + //)) + 3L (L + 1)] (4) 

For fixed A, [i one gets the spectrum of a rotational band, but going as L (L + 1) rather than I (I + 1). 

We thus see that the Q ■ Q interaction has been used to make more transparent how nuclear collectivity arises 
despite the fact that the shell model with its implications of single particle motion seems to work. It should be added 
that one can regard Q ■ Q as the long range part of a more realistic interaction. This has been shown by many authors 
and we recommend Talmi's book [pf for a nice detailed derivation of this. 



II. MOMENTUM DEPENDENT Q.Q INTERACTIONS AND THE STABILITY PARADOX 



It is of general interest to know how momentum dependent interactions affect nuclear deformations. For example 
the Skyrme interaction contains the momentum dependent terms k' 2 S (ri — + S (ri — k 2 and k'-<5 (ri — r^) k 
p5| |l7| . The presence of those terms leads to an effective mass m* < m in the nucleus. The rationale given to 
those momentum dependent terms is that they simulate the finite range nature of the nucleon-nucleon interaction. 
However, since the Skyrme interactions are phenomenological they may be simulating more fundamental momentum 
dependence of the nucleon-nucleon interaction. 

The effect of the momentum dependent terms is known for the energy of the isoscalar giant quadrupole resonance 



E* 



V2huj, — 
V m* 



(5) 



The effect on the ground state deformations is not so well known. 

Another example of momentum dependence is the Hamiltonian used by Elliott to obtain SU(3) results 
Hamiltonian is symmetric in position and momentum coordinates 
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Q r + Q p 



Q r + Q p 



(6) 



Ql = r 2 Y 2 ^ (O r ) ; Ql - b 4 p 2 Y 2 ^ (Q p ) ; (6 2 = h/mu) . (7) 
We also define for convenience 

= \ (Q r + Q p ) ■ (8) 



^E 



The momentum terms prevent AiV = 2 mixing between major shells. Although the original intention for this 
Hamiltonian is for valence nuclcons in a given major shell, let us use it here as an extreme example of a momentum 
dependent Hamiltonian which is symmetric in r and p and perform a Hartree-Fock calculation with it. 

In general, for any interaction in a deformed Hartree-Fock calculation the condition for stability can be formulated 
as follows Ijl 911 

^E(^))= ' (9) 

where Q p is the quadrupole moment operator in momentum space 



2 



However, since the Hamiltonian is symmetric in r and p, so is the Wigner distribution function [|20|. This leads to 
the conclusion (J^. Q r ) — 0, i.e., the expectation value of the usual quadrupole moment also vanishes. But then this 
goes against the common belief that the Elliott model supports rotational motion. 

A way out of this dilemma is to do, as we have done, go back to the r space interaction Q r ■ Q r and note that we 
can still get the SU(3) results by including in the single-particle energies, not only the contribution from the i = j 
term in the sum 

\Y.Q r {i)-Q r (j) 

i,j 

but also the particle-core interaction |2^j2^] . In these references it was noted that Q r ■ Q r and Q E ■ Q E do not give the 
same results in a valence space even if one uses harmonic oscillator wave functions unless the particle-core interaction 
is included. We can write 

5£qW-Q(j) = £Q(0-Q0") + 5£qW-Q« • (ii) 

ij i<j i 

For the first term one does get the same result for Q = Q E and Q — Q r in a given N shell. However this is not the 
case for the second i = j term. Using Q = Q r one only gets 2/3 of the value obtained with Q = Q E . We have however 
shown ^l],^2| that the remaining 1/3 comes from the (exchange) interaction of the valence nucleon with the core. 

Just to clarify things, the diagonal part of Q r ■ Q r and the particle-core interaction are the terms that give the 
single particle splittings Anu,£') for different £ values in a major shell. In order to get the SU(3) results of Elliott 
one must have such a single particle splitting as well as the two-body Q ■ Q interaction between the valence nucleons. 
The main difference then between Q E ■ Q E and Q r ■ Q r is that with the former all the single particle splitting comes 
from the i = j part of Q E ■ Q E , whereas with Q r ■ Q r two thirds comes from the i = j part and one third from the 
particle-core interaction. In more detail the expressions for the single particle splitting between N£ and N£' states 
are || 

A 1 N{i>e) =2x[£(£ + l)-£'(£' + l)} , 

A 2 W) =x[^+l)-^' + l)] , 
where the AL,« comes from the i = j part of Q r ■ Q r and A 2 N / e en comes from the particle-core interaction. 



III. 1(1+1) VS L(L+1) FOR ROTATIONAL BANDS 

In this Section we stay in the valence space and make a comparison of Elliott's SU(3) results f|-|(| with results of 
the rotational model |ll|. We feel that it has not been sufficiently emphasized that one finds cases where these two 
models yield different results for the behavior of rotational bands. 

In the rotational model the formula for the energy of a state in a rotational band with total angular momentum / 
is given by [pi 



h 2 

Ei = E a H 

21 



1(1+1) + Sk,i/2 (-l/ +1/2 (/ + 1/2) a] (12) 



where a is the decoupling parameter given by a = — (K =1/2 | J+| K = 1/2) and where if \K) = J2j ^j,K4>j,K then 

\K) = Y lj C j ,K{-l) j+K <t> j ,-K. 

As noted by J.P. Davidson |2^] for an odd-odd nucleus with K n = ±1/2, K p — +1/2, (K — 0), there is an 

additional term (— 1) +1 a p a n SK,o ( < 5r" p ,i/2 < 5r"„.-i/2 + &K p ,-i/2&K n ,i/2) (See also a recent review by Jain et al. pjf). 
For even-even nuclei, and for odd-even and even-odd nuclei with K + 1 1/2, one gets the familiar I (I + 1) spectrum 

& 

It is generally thought that the Elliott SU(3) model also gives an 1(1+1) spectrum. A careful reading of the 
papers however shows that one really gets an L (L + 1) spectrum, where L is the orbital angular momentum |}]-||. 
In this work we wish to point out that to fully convey the similarities and differences of the Elliott model and the 
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rotational model, one should consider not only even-even nuclei, but also even-odd (odd-even) and especially odd-odd 
nuclei. The latter are usually not considered in the standard textbooks. 

The point we make in this section is that the spectrum in the SU(3) model is indeed L (L + 1) in all of the above 
instances. This leads to several cases: 

(a) Rotational bands where the spin S is equal to zero. In this case I = L and one gets an 1(1+1) 
spectrum. 

(b) K = 1/2 bands of even-odd (odd-even) nuclei. Here again the SU(3) spectrum is of L (L + 1) type. The 
Elliott formula is the same as the rotational formula |l|] including the decoupling term evaluated in the asymptotic 
limit. Again, we have consistency between the two approaches. 

(c) Odd-odd nuclei. Here again the SU(3) spectrum is L (L + 1). There are some cases here where one gets a 
behavior which is not consistent with the rotational formula. 

We have performed shell model calculations with all possible configurations in a given major shell using the inter- 
action J2i<j QW ' Q{j) where, in order to get Elliott's SU(3) results we must, as mentioned in Section II, also add 
single-particle splittings. Therefore, in the Is — Od shell we have eod — ei s = 18% and in the lp — 0/ shell we have 
£o/ — tip = 30x , where \ — (56 4 /327r) \ with b the oscillator length parameter b 2 = h/mco [ pip2l . 

The same results can of course be obtained from the Elliott SU(3) formula for the energies given in Eq.(|j) 

One has the further rules |^J: 

Let A be the maximum of A and /i, and jl the minimum. Then Kl — fl, p. — 2, 1 or and 

• L = K L ,K L + 1, ...,K L + A when K L £ 0; 

• L = A, A - 2, 1 or when K L = 0. 



A. A brief look at K = 1/2 bands 

Let us be specific and discuss 19 F and 43 Sc. We consider in each case three valence nucleons beyond a closed shell. 
In 19 F the particles are in the Is — Od shell, whereas in 43 Sc they are in the lp— 0/ shell. The energy levels of the lowest 
bands are given in Table I for the two cases. The results for the two nuclei are striking but different. In 19 F, the lowest 
state is a / = 1/2+ singlet, and at higher energies we get degenerate pairs (3/2 + , 5/2 + ), (7/2+, 9/2 + ), (11/2+, 13/2+). 
In 43 Sc the ground state is degenerate, and the degenerate pairs are (l/2~,3/2~), (5/2 - ,7/2~), (17/2~, 19/2~). 

If we look at the rotational formula, we find that these results are consistent with a decoupling parameter a = +1 
for 19 F and a = — 1 for 43 Sc. It is easy to show that these are precisely the a— values one obtains with asymptotic 
Nilsson wave functions. In both cases the odd particle will be in a A = 0, E = 1/2 state in the asymptotic limit. From 

the definition of K, the state 



A = E = 1/2) can be shown to be equal to |A = E = -1/2) , where (— l) w 

is (+1) for an even-parity major shell and (-1) for an odd-parity one. Hence: 

a = (-If (E = +1/2 | J+| E = -1/2) = {-If (13) 

It has long ago been noted by Bohr and Mottelson |l3| that a = +1 corresponds to weak coupling of the odd 
particle to / = + , 2 + , 4 + , ... states, whereas a — — 1 corresponds to weak coupling to / = 1, 3, 5 states. In the context 
of the SU(3) model, we would say that the 19 F states have an L(L + 1) spectrum with only even L's allowed, and 
that 43 Sc has an L(L + 1) spectrum with only odd L's allowed. It should be emphasized that the purpose of Table 
I is to compare the results of the Q ■ Q interaction with the rotational model. To compare with experiment, more 
realistic interactions including in particular spin-orbit may be required. In this respect, the results in Table I can be 
considered as the asymptotic limit of large deformation. 

At any rate, we have shown that the Q ■ Q interaction gives the same results for these two K = 1/2 bands as does 
the rotational formula with asymptotic Nilsson wave functions. 

Although it is not our intention in this work to fit experiment, rather we wish to work out the consequences of our 
model, we cannot resist discussing an interesting result for 19 F. There is a beautiful example of weak coupling in 19 F 
corresponding to a = +1, but it is not the positive parity band discussed above but rather a negative parity band 
corresponding to a 4-particle 1-hole configuration. This has been discussed by Nazarewicz et al. [p5[ . In the asymptotic 
limit the proton hole would be in the Nilsson orbit [101] \ which would yield a = 0. However, the dominance of the 
P1/2 due to spin-orbit interaction causes a to be very close to one. The spectrum looks like the weak coupling of a 
P1/2 hole to the states of the ground state rotational band of 20 Ne. 
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B. Odd-odd nuclei, e.g. 22 Na 

In Table II we show a fairly detailed list of energy levels for the odd-odd nucleus 22 Na obtained with the Q ■ Q 
interaction. We show T = and T = 1 states in separate columns. We use the same parameters as in 19 F just to 
bring out some similarities. If one is interested in a best fit, one should of course have an A dependence in \- 

An striking feature in this table is that many states with diverse VT assignments are degenerate. This is due to 
the prevailing SU(4) ® SU(3) symmetry. We shall come back to this point at the end of this section but first we 
compare with the rotational model. 

We have underlined T — and T = 1 rotational bands, and will now discuss them in more detail. Note that the 
ground state consists of two degenerate states, one with / = 1 + T = and the other with I = + T = 1. Both states 
have L — and the simple spin-independent interaction gives the same energy for S = and S = 1 . Let us first look 
at the underlined T = 1 states. The ground state is I = + , the 2 + state is at 1.588, the 4 + is at 5.293, etc. If we 
follow the rotational sequence I = + , 2 + , 4 + , ... we see a simple rotational behavior: 



E{I) - E (0+) = AI (I + 1) ; with A = 

There is nothing new here. 



h 2 _ E{2+)-E(Q+) 



21 6 

The excitation energies of the underlined T = 0, I = 1 + , 2 + , 3 + , 10 + states follow the sequence 



E(2+)-E(l+) 

E*(I) = E(I) - E(lt) = AI(I + 1); with A = — - 1 1 ' 
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However the rotational model formula (Eq. (|l2|)) would give a different energy level spacing 

E*(I) = E{I) - E(l+) = A' [1(1 + 1) - 2] ; with A' = - ^ ~ - ^ \ 

which is not followed by the above mentioned energy levels in the table. Thus, for the case of T = states in odd-odd 
nuclei we get a clear difference between the rotational formula and SU(3). 

We gain further insight by examining the degeneracies associated with the T = underlined states in Table II, i.e., 
those with energy AI(I + 1). The even / states up to I — 8 are doubly degenerate whereas the I — 10 and the odd-/ 
states are singlets. This suggests that there are two bands for which the states with the same I values are degenerate. 
One band is a K = 2 band with all values of I from 2 to 10, and there is nothing anomalous about it. The other band 
consists of states of angular momentum 1,2,4,6 and 8. For the latter band, the orbital angular momentum of the states 
are 0,2,4,6 and 8 respectively, and they all have 5=1. Their energies can be fit to the formula E*(I) — AL(L + 1) 
rather than AI(I + 1), such that only even L contribute. 

Let us now discuss on the basis of SU(3) the degeneracy observed in Table II, limiting the discussion to T = and 
T = 1 states (similar arguments would follow for higher isospin states). For 22 Na we have 3 protons and 3 neutrons 
in the N — 2 shell. The ground state will correspond to maximal spatial symmetry in coordinate space, i.e., to the 
(A, (i) — (8, 2) representation of the [/] = [4, 2] partition. The possible Kl, L values are 

• K L =0 : L = 0,2,4,6,8 

• K L = 2 : L = 2, 3, 4,5, ...,10 

Since the energy does not depend on Kl (see Eq.([|)), states with equal L— values and different Kl are degenerate. 
In addition, antisymmetry in spin-isospin space demands that S = when T = 1, and S = 1 when T = 0. Therefore, 
T = 1 states with I = L, and T = states with I = L, L + 1,L — 1 will be degenerate. Thus for the ground state 
[4, 2] (8, 2) representation we find with the same excitation energy E* [L] = i\L (L + 1) , the following states: 

• For L = 2, 4, 6, 8 

— two states, each with T = and I — L, L ± 1, and 

— two states, each with T = 1 and I = L. 

• For L = 3,5, 7, 9, 10 

— one state with T — and I = L, L ± 1, and 
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— one state with T = 1 and I = L. 

This explains why the state with excitation energy 1.588 MeV (i.e., L = 2) appears twice for T = and I = 1, 2, 3 
and twice for T = 1 and 7 = 2. The state with excitation energy 3.176 MeV (i.e., L = 3) appears once with 
T = 0, 1 = 2, 3, 4 and once with T = 1, 1 = 3. The next state of this representation is at 5.293 MeV (i.e., L = 4), and 
according to the above discussion should appear twice for T = and I = 3, 4, 5 and twice for T = 1 and 7 = 4. 

Surprisingly though there are more degeneracies at 5.293 MeV, namely one T = state with 7 = 3 and three 
T = 1 states with 7 = 2,3, 4. This is due to an additional degeneracy of [4, 2] (8, 2) L = 4 with the [4, 1, 1] (9, 0) L = 3 
state that corresponds to total symmetry in spin-isospin space (i.e., T = 0, S = or T = 1,S = 1). The allowed 
Kl,L values in this representation are Kl = 0, L = 1,3,5,7,9, and the excitation energy for a given L— value, 
E* [{§Q)L]~E* GS [0] = x [24 + 3L [L + 1)] , is shared by states with T = 0, 1 = L and with T = 1,1 = L, L±l. Thus the 
states that belong to this representation for L = 1 are the states at 2.647 MeV with T = I = 1 and T = 1 1 = 0, 1, 2, 
while for L = 3 we get the states at 5.294 MeV previously discussed with T = I = 3 and with T = 1 7 = 4, plus 
the states 7 1 = 1 7 = 2, 3. The next states with L = 5 (75* = 10.055) appear for T = I = 5, T = 1 I = 4, 5, 6, 
and are degenerate with lower angular momentum states corresponding to Kl = 1, L = 2 in the (6, 3) representation 
(The L = 1 states in this last representation have 9 MeV excitation energy). By the same token one can explain the 
remaining energy levels in the table. To the best of our knowledge there has been no discussion previously of these 
additional degeneracies corresponding to different values in different representations. 

Experimentally the lowest band in 22 Na has K = 3. In the Nilsson scheme when we fill the lowest N = 2 K = 1/2 
level we reach 20 Ne. To form 22 Na we put the odd neutron and odd proton into the K = 3/2 level and the Gallagher 
Moszkowski rule favors K = 3 over K = for the lowest band p(|. However, with pure Q ■ Q as seen in Table II, the 
lowest 7 = 3 + state comes out at 1.588 MeV of excitation. Clearly, the limitations of the SU(3) formula for energy 
levels are more apparent when one examines the spectra of odd-even and odd-odd nuclei rather than limiting oneself 
to even-even nuclei. The main deficiency is the absence of a spin-orbit interaction in the SU(3) model. For even-even 
nuclei the lowest lying levels are dominantly S = states but this is not the case for T = states of odd-odd nuclei. 



IV. CONVERGENCE WITH A Q.Q INTERACTION 



We now come to the last topic of this work - how to achieve convergence with a Q • Q interaction. When we remove 
the momentum dependent terms from Q ■ Q and go back to the r space Q r ■ Q r , we can mix major shells. Therefore, 
this coordinate space interaction combines the advantages of the algebraic Elliott's model, with the possibility of 
studying intruder states and admixtures of higher shells in the ground state. 

There is a general belief that if we go to larger and larger spaces the results will diverge. This is because the Q r ■ Q r 
interaction increases its strength as the distance between two nucleons increases. However, we wish to show in this 
work that we can get around this problem. Just as people do all the time in G-matrix calculations, we can modify 
the interaction as we change the model space. If we allow up to nhu excitations in the model space we can make the 
strength of the Q r ■ Q r depend on n 

V = -^yQ r ■ Q r . 

A simple scheme for getting convergence is to demand that the energy of the I = 2\ state in an open shell nucleus 
comes out correctly (i.e., agrees with experiment) for each n. We have applied this scheme to the 7 = + , 1 + and 2 + 
(all T = 0) states in 8 Be. These are shown in Tables III, IV, and V, respectively. In each table the results are given 
in a OTilu space (n = 0), (0 + 2)hu> space (n = 2), and (0 + 2 + A)%uj space (n = 4). In order to get the 7 = 2^ state 
to come out correctly we had to steadily decrease x with increasing n. The values for n = 0, 2, and 4 are respectively 
0.5216, 0.4119, and 0.3369 MeV fm^ 4 . 

The calculations have been performed using the OXBASH program p7| . In this program the effects of spurious 
states are removed using the Gloeckner-Lawson method [p8| . In this method the spurious states are pushed up to a 
very high energy. This is done by diagonalizing the modified Hamiltonian for a system of A nucleons 



H' = H SM + A 



CM 



1 3 

— TO^4tt> 2 R 2 hu 



2 Am 2 2 



(14) 



where Hsm is the usual shell model Hamiltonian and P and R are the position vectors for the center of mass. By 
making Xcm very large, the spurious states corresponding to center of mass motion will be pushed up to a very high 
energy. They will be well separated from the lower lying physical states. 
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Let us first discuss the calculations up to 2hui, the results of which are given in the tables III, IV and V. These 
excitations are of two types, we can either excite one particle through 2 major shells (lp — lh) or we can excite two 
nucleons through one major shell (2p — 2h). The (lp — lh) excitations can be regarded as giant resonance states built 
upon the valence states as have indeed been studied in the context of symplectic symmetry . The giant resonance 
states are generally at high excitation energies but their admixtures into the valence states can lead to important 
effects. For example, the isoscalar E2 effective charge needed to fit experiments in valence space calculations is 
about a factor two larger than the bare charge and this can be understood in terms of small admixtures of the 2Hlo 
(lp — Instates into the basis valence states. 

We next consider the 2p — 2h states. In some nuclei some states of this type come much lower in energy than the 
2hu) estimate. They may intertwine with low lying valence states. In such cases these 2p— 2h states or more generally, 
np — nh states, are called intruder states. In previous calculations p9| we have shown that there are such low lying 
intruders in 10 Be, 12 C and 16 0, but actually not in 8 Be. 

The reason that in some nuclei the intruders come low in energy and in others they do not, can be understood 
in the context of the Nilsson model. The optimum way to have a low lying intruder is to lift two nucleons from a 
Nilsson orbit whose energy increases with deformation (upgoing level) and to put them into a Nilsson orbit whose 
energy decreases with deformation (downgoing level). For 10 Be, 12 C and 16 we can remove nucleons from upgoing 
levels in the Op— shell, but in the case of 8 Be the valence levels are depleted and we would have to remove nucleons 
from a downgoing level. 

Before discussing the energy levels we will make some comments about the occupancies. It had been previously 
noted by Fayache et al. J3(| that if one allows only 2hui admixtures with the Q ■ Q interaction then the shell model 
matrix reduces into two parts. First, we have states which are admixtures of valence states (no particles excited to 
higher shells) and states in which one nucleon is excited through two major shells, and second we have states in which 
two particles are excited through one major shell. There is no mixing in between the (Op — Oh + lp — lh) states and 
the 2p — 2h states. This is because of the 'parity' selection rule that Q ■ Q cannot excite two nucleons from an even 
(odd) parity major shell to an odd (even) major shell. 

If we look at Tables III, IV, and V under the columns (0 + 2) hui we see several states which have 100% 2hcu 
configurations. These are precisely the 2p — 2h states which are uncoupled from the rest. When we allow 4hui 
admixtures the 2p — 2h states admix with states in which one nucleon is excited through one major shell and the 
other through three major shells. We therefore no longer get 100% 2hui for these states, but they are still uncoupled 
from states which have mainly valence configurations. 

In more detail, for / = + in Table III, in the (0 + 2) hui column the ground state has 21.77% 2hui admixtures. 
The next four states have smaller 2hui admixtures, but the state at 30.951 MeV has an 88.89% admixture. This is 
clearly a state which is dominantly a lp — lh state in which a nucleon has been excited through 2 major shells. This 
state can be easily distinguished from the 2p - 2h states at 32.214 MeV and 34.409 MeV. The latter have 100% 2hcu 
configurations. As mentioned before there is no mixing, at the 2Hlo level between the (Op — Oh + lp — lh) states and 
the 2p — 2h states via a Q ■ Q interaction. 

States with 100% occupancy are also seen in Tables IV and V corresponding to / = 1 + and 2 + . These all lie in the 
energy range 30-35 MeV, and are clearly 2p— 2h states. As mentioned before there are no low lying intruders in 8 Be. 
These 2p — 2h states are close in energy with the (dominantly) lp — lh states for 8 Be. 

When we go to the (0 + 2 + 4) hui column the results are not so striking - we do not see 100% of any configuration. 
But we can still distinguish the (Op — Oh + lp — lh) states from the (2p — 2h + 2p — 2h <£> lp — lh) states. For the 
latter the sum of the 2hui + Ahui occupancy should be 100%. For the I — + states in Table III this is not the case 
for the state at 26.225 MeV (71.5%) but it is the case for states at 32.085 MeV and 34.200 MeV. Similar states can 
be found for I = 1 + and 2 + , i.e. states where the 2hui + Afrui occupancy is 100%. This confirms that with Q ■ Q one 
gets a separation of different classes of states. 

We now discuss the energy levels as a function of n, the number of hco excitations. While there are rather large 
deviations in going from n = to n = 2, there is excellent convergence for many states in comparing n = 2 and 
n = 4. The percentage deviation for the I — + case for the first 4 states (up to about 20 MeV) are respectively 
0.77, 2.22, 2.75, and 3.93%. For the next state, which is at 26.225 MeV in the (0 + 2 + 4)hu> calculation, the percentage 
deviation is large 18.02%. But if we examine this state we see that it is mainly a 2hui excitation state, not present 
in the Ohui calculation. It corresponds to a lp — lh excitation through 2 major shells, as opposed to many other 2Hlo 
states which correspond to two nucleons being excited through one major shell. 

We see that the convergence between (0 + 2) hui and (0 + 2 + 4) hui holds to surprisingly high energies. For example 
for / — 1 + at the (0 + 2) hui energy level there are five consecutive 2p — 2h states between 30.722 MeV and 34.503 
MeV. The % deviations in the energies relative to the (0 + 2 + 4) hw calculations are respectively 10.13, 9.19, 4.84, 
2.01 and 1.87%. 
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Table I. Energy levels (MeV) of excited states corresponding to the K = 1/2 ground state bands in 19 F and 43 Sc 
with the — xQ • Q interaction. 



19p a 43 Sc b 



r 




r 


E* 


(1/2)+ 





(1/2)" 





(3/2) + 


1.588 


(3/2) 





(5/2) + 


1.588 


(5/2)" 


0.679 


(7/2)+ 


5.295 


(7/2)" 


0.679 


(9/2)+ 


5.295 


(9/2)" 


1.900 


(11/2) + 


11.118 


(11/2)" 


1.900 


(13/2) + 


11.118 


(13/2)" 


3.664 






(15/2)" 


3.664 






(17/2)" 


5.971 






(19/2)" 


5.971 


a p Qr 19p we uge . 


= 0.1841 (x = 0.0882) 






b For 43 Sc we use \ 


= 0.0294 {x = 0.0218) 
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Table II. T = and T = 1 energy levels (MeV) of 22 Na calculated with the — xQ ■ Q interaction a . Only the first six 
levels for each I* are shown. 



T = states 



T = 1 states 



0+ 



8.999 
12.176 
12.176 
13.235 
16.410 
16.411 



0.000 
2.647 
8.999 
9.000 
12.176 
12.176 



T = states 



T = 1 states 



6+ 



7.941 

11.117 

11.117 

14.824 

16.411 

16.411 



10.059 
11.117 
11.118 
16.412 
16.412 
16.412 



3+ 



5+ 



0.000 


2.647 


7+ 


11.117 


14.823 


1.588 


8.999 




11.117 


16.940 


1.588 


8.999 




14.823 


19.587 


2.647 


10.059 




16.941 


19.587 


9.000 


10.059 




19.058 


19.587 


9.000 


10.059 




19.059 


19.587 


1.588 


1.588 


8+ 


14.823 


16.941 


1.588 


1.588 




19.058 


19.058 


3.176 


2.647 




19.059 


19.059 


8.999 


5.294 




22.763 


22.767 


10.059 


9.000 




23.292 


23.293 


10.059 


9.000 




23.293 


23.293 


1.588 


3.176 


9+ 


19.058 


23.822 


1.588 


5.293 




19.058 


25.939 


3.177 


10.058 




23.822 


26.470 


5.294 


10.058 




25.940 


27.527 


5.294 


11.646 




26.469 


27.527 


5.294 


11.646 




27.528 


29.644 


3.176 


5.293 


10+ 


23.823 


25.942 


5.293 


5.293 




29.117 


29.117 


5.293 


5.293 




30.705 


30.706 


7.941 


10.059 




32.293 


32.294 


11.647 


11.647 




33.881 


32.294 


11.647 


11.647 






33.882 


5.294 


7.941 








5.294 


10.059 








7.941 


13.763 








10.059 


13.763 








11.118 


13.763 








11.118 


13.763 









a In this table and in the following tables, the same value of \ ( an d of \) was used for 22 Na as for 19 F. 
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Table III. The energies E* (MeV) and percentages of 2hui and 4hui occupancies of the J = + T = states in 8 Be in 
three different model spaces. The last column gives the % energy deviation between the (0 + 2) %lo and (0 + 2 + 4) hui 
calculations. 



Ohuj (0 + 2) Tiuj (0 + 2 + 4) tuj % deviation 

X = 0.5216 x = 0.4119 x = 0.3369 

E* E* %2huj E* %2Tioj %ATioj 



0.00 


0.000 


21.77 


0.000 


24.60 


10.74 




9.12 


11.458 


10.28 


11.370 


10.90 


5.69 


0.77 


12.16 


16.231 


1.91 


15.879 


2.07 


3.55 


2.22 


15.20 


18.352 


2.96 


17.860 


2.53 


3.16 


2.75 


17.23 


20.147 


2.62 


19.385 


2.10 


2.98 


3.93 




30.951 


88.89 


26.225 


50.93 


20.56 


18.02 




32.214 


100.00 


29.701 


77.30 


19.41 


8.46 




32.396 


92.36 


32.085 


86.13 


13.87 


0.97 




34.409 


100.00 


34.200 


86.82 


13.18 


0.61 




38.702 


91.66 


35.928 


70.75 


15.45 


7.72 
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Table IV. Same as Table III but for the J = 1+ states. 



Hlo (0 + 2) hco (0 + 2 + 4) huj % deviation 



E* 


E* 


%2hw 


E* 


%2hco 


%4hio 




9.12 


11.458 


10.28 


11.371 


10.90 


5.69 


0.77 


11.15 


13.734 


8.85 


13.587 


8.92 


4.87 


1.08 


15.20 


18.353 


2.96 


17.861 


2.53 


3.16 


2.75 


17.23 


20.148 


2.62 


19.385 


2.10 


2.98 


3.94 




30.722 


100.00 


27.896 


78.69 


21.31 


10.13 




32.992 


100.00 


30.215 


79.30 


20.70 


9.19 




QQ on A 


1UU.UU 


oL. 100 


OZ.00 


1 ( .40 


4.84 




34.409 


100.00 


33.731 


86.28 


13.72 


2.01 




34.503 


100.00 


33.869 


83.31 


16.59 


1.87 




35.357 


99.84 


34.200 


86.81 


13.19 


3.38 




37.666 


92.98 


35.707 


74.99 


14.33 


5.49 




38.702 


91.65 


35.927 


70.75 


15.44 


7.72 




40.174 


100.00 


38.724 


85.12 


14.80 


3.74 




40.485 


99.97 


39.313 


86.37 


13.63 


2.98 




40.706 


100.00 


39.331 


86.38 


13.63 


3.50 
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Table V. Same as Table III but for the J = 2+ states. 



Hlo (0 + 2) hco (0 + 2 + 4) huj % deviation 



E* 


E* 


%2hio 


E* 


%2hio 


%4hio 




3.04 


3.042 


21.42 


3.038 


23.81 


9.94 


0.13 


9.12 


11.458 


10.28 


11.370 


10.90 


5.69 


0.77 


11.15 


13.734 


8.85 


13.586 


8.92 


4.87 


1.09 


12.15 


16.197 


8.99 


15.880 


2.07 


3.55 


2.00 


14.18 


16.231 


1.91 


15.949 


8.31 


4.21 


1.77 


1 £ on 


lo.oOo 


z.yo 


1 1 S^O 

1 ( .ooy 


Z.OO 


o.lu 


O 11 


17.23 


20.148 


2.62 


19.385 


2.10 


2.98 


3.94 




31.711 


82.81 


27.146 


51.37 


20.15 


16.82 




32.214 


100.00 


30.215 


79.30 


20.70 


6.62 




32.992 


100.00 


31.714 


80.12 


18.91 


4.03 




33.962 


99.68 


32.085 


86.13 


13.87 


5.85 




34.410 


100.00 


33.869 


83.31 


16.60 


1.60 




35.357 


99.84 


34.201 


86.81 


13.19 


3.38 




37.668 


92.98 


35.707 


75.00 


14.33 


5.49 
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